JUN  79  H I SANCER#  A 0 VARVATSIS*  S SIEGEL  F29601-78-C-0098 


UNCLASSIFIED  AFWL  -TR-78-172  NL 


DDC  file  copy  ft  A 071 42 


AFWL-TR-7  8-172 


AFWL-TR 

78-172 


PSEUDOSYMMETRIC  EIGENMODE  EXPANSION  FOR 
THE  MAGNETIC  FIELD  INTEGRAL  EQUATION 


AND  SEM  CONSEQUENCES 

Maurice  I.  Sancer 
A.  D.  Varvatsis 
Scott  Siegel 

TDR,  Inc. 

Marina  del  Rey,  CA  90291 


June  1979 


Final  Report 


Approved  for  public  release;  distribution  unlimited. 

This  research  was  sponsored  by  the  the  Defense  Nuclear  Agency 
under  Subtask  R99QNXEB200 , Work  Unit  71,  Work  Unit  Title: 
Hardness  Surveillance 


Prepared  for 
Director 

DEFENSE  NUCLEAR  AGENCY 
Washington,  DC  20305 

AIR  FORCE  WEAPONS  LABORATORY 
Air  Force  Systems  Command 
Kirtland  Air  Force  Base,  NM  87117 


DDC 


! 


*»  07  17  oil 


AFWL-TR-78-172 


4 


This  final  report  was  prepared  by  TDR  Incorporated,  Marina  del  Rey, 

California,  under  Contract  F29601-78-C-0098,  Job  Order  WDNF3206  with  the  Air 
Force  Weapons  Laboratory,  Klrtland  Air  Force  Base,  New  Mexico.  Captain  Howard 
G.  Hudson  (ELTI)  was  the  Laboratory  Project  Officer-1 n-Charge. 

When  US  Government  drawings,  specifications,  or  other  data  are  used  for  any 
purpose  other  than  a definitely  related  Government  procurement  operation,  the 
Government  thereby  Incurs  no  responsibility  nor  any  obligation  whatsoever,  and 
the  fact  that  the  Government  may  have  formulated,  furnished,  or  In  any  way 
supplied  the  said  drawings,  specifications,  or  other  data.  Is  not  to  be  regarded 
by  Implication  or  otherwise,  as  In  any  manner  licensing  the  holder  or  any  other 
person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture,  use, 
or  sell  any  patented  Invention  that  may  In  any  way  be  related  thereto. 

This  report  has  been  authored  by  a contractor  of  the  United  States  Government. 
Accordingly,  the  United  States  Government  retains  a nonexclusive,  royalty-free 
license  to  publish  or  reproduce  the  material  contained  herein,  or  allow  others 
to  do  so,  for  the  United  States  Government  purposes. 

This  report  has  been  reviewed  by  the  Office  of  Information  (01)  and  Is 
releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS,  It 
will  be  available  to  the  general  public.  Including  foreign  nations. 

This  technical  report  has  been  reviewed  and  Is  approved  for  publication. 


HOWARD  G.  HUDSON 
Captain,  USAF 
Project  Officer 


hlef.  Technology  Branch 


FOR  THE  COMMANDER 

/fitiyvci-i)  & /Q&urQ, 

DONALD  A.  DOWLER 
Colonel , USAF 

Chief,  Electromagnetics  Division 


DO  NOT  RETURN  THIS  COPY.  RETAIN  OR  DESTROY. 


o " 


Vb 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Enterad) 


REPORT  DOCUMENTATION  PAGE 


)i/»C  READ  INSTRUCTIONS 

AUC- BEFORE  completing  form 

2.  GOVT  ACCESSION  NO.  3.  RECIPIENT'S  CATALOG  NUMBER 


AIW-rR-78-172 


U.  TITLE  fand  Sub 


5EUD0SYMMETRIC  EIGENMODE  EXPANSION  FOR  THE  / 
^GNETIC  FIELD  INTEGRAL  EQUATION  AND  SEM  ^ 
JNSEQUENCtS#  ^ ^ ==•  ( 


S.  TYPE  OF  REPORT  ft  PERIOD  COVERED 


17.  AUTHORfa) 


(k\  Maurice  I./Sancer,  A.  D./Varvatsls 
KH  Scott^Slegel  i- 1 *m 


Final 


8.  contract  or  grant  number^*; 

^296^1  -78-C-py98 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  ft  WORK  UNIT  LUMBERS 


9.  PERFORMING  ORGANIZATION  NAME  AND  AOORESS  10.  PROGRAM  ELEMENT.  PROJECT,  ' 

. AREA  ft  WORK  UNIT  LUMBERS 

Inc-  62704H  - DNA 

P.  0.  Box  9695  64747F  - AFWL 

Morlwa  de4Hteyv  CA  90291  Le>s<q-^;f*dL  if/)  HDNE3206 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  ^ ~S.  RCRO'RT'PRTB  

Director  " x t / / ( J ] JunJ-4979  f 

Defense  Nuclear  Agency  * - y " ^ -if- number  bp  phubi  

Washington,  DC  20305  58 

U.  MONITORING  AGENCY  NAME  ft  ADORESSflf  dilt'rtnt  trom  Controlling  Office)  15.  SECURITY  CLASS,  (of  this  report) 

Air  Force  Weapons  Laboratory  (ELTI)  min  jccrcrcn 

Klrtland  Air  Force  Base,  NM  87117  UNLLAbMritu 


4*!  (j  i 


NtPUM  UNIL- 

11)  ihinf-4379 


UNCLASSIFIED 


15a.  DECLASSI  FI  CATION/ DOWN  GRADING 
SCHEDULE 


I 16.  DISTRIBUTION  STATEMENT  (ol  thin  Report) 


Approved  for  public  release;  distribution  unlimited. 


p- 


17.  DISTRIBUTION  STATEMENT  (ot  the  abstract  entered  In  Block  20,  II  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


This^research  was  sponsored  by  the  Defense  Nuclear  Agency  under  Subtask 
R99QI|fXEB200,  Work  Unit  71,  Work  Unit  Title:  Hardness  Surveillance 


1 19.  KEY  WORDS  ( Continue  on  reveres  aide  It  necessary  and  Identify  by  block  number) 


Electromagnetic  Scattering 
Perfectly  Conducting  Bodies 
Perfectly  Conducting  Sphere 
Magnetic  Field  Integral  Equation 


Eigenmode  Expansion  Method 
Singularity  Expansion  Method 
SEM  Patch  Zoning 


20.  ABSTRACT  ( Continue  on  reverse  aide  It  necessary  and  Identify  by  block  number) 

An  eigenmode  expansion  for  the  magnetic  field  Integral  equation  (MFIE)  Is 
derived  which  eliminates  the  requirement  that  an  adjoint  solution  be  explicitly 
sought.  Instead,  an  orthogonality  relation  Is  derived  which  only  Involves  the 
elgenmodes  of  the  original  MFIE  operator.  To  promote  confidence  In  the 
validity  of  the  resulting  expansion,  two  analyses  based  on  this  expansion  are 
presented  which  lead  to  known  results.  First,  the  expansion  Is  applied  to  the 
problem  of  determining  the  surface  current  density  Induced  on  a perfectly 
(over) 


DD  | JAN *7*  1473  EOITION  OF  I NOV  «S  It  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  flW>»n  Dmtm  Enfrtd) 


M mb 


Sr 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfWimn  Data  Entered) 


20.  ABSTRACT  (Cont'd) 

conducting  sphere  by  a plane  wave  and  the  known  solution  for  this  problem  Is 
duplicated  by  the  expansion.  The  second  analysis  shows  that  for  a general 
perfectly  conducting  body,  the  eigenmode  expansion  coefficient  numerator 
evaluated  at  the  purely  Imaginary  frequency  corresponding  to  an  Interior  reso- 
nance Is  zero.  This  result  Is  necessary  In  order  to  relate  the  eigenmode 
expansion  to  SEM.^^ 

Viewing  the  SEM  as  a change  of  representation  of  the  eigenmode  expansion 
Intended  to  facilitate  the  Inverse  Laplace  transform  of  that  expansion  and 
taking  advantage  of  our  detailed  sphere  calculations,  we  obtained  an  Important 
SEM  result.  We  found  that  MFIE  class  2 coupling  coefficients  give  the  wrong 
answer  for  the  sphere.  This  result  caused  us  to  examine  class  2 coupling 
coefficients  corresponding  to  the  electric  field  Integro-dlfferentlal  equation 
(EFIDE).  We  examined  the  symmetric  eigenmode  expansion  corresponding  to  that 
equation  and  found  that  EFIDE  class  2 coupling  coefficients  also  give  the  wrong 
answer  for  the  sphere.  This  detailed  examination  of  the  sphere  solution  led  us 
to  postulate  a set  of  SEM  assumptions  that  have  potential  application  to  a 
general  class  of  closed  surfaces  and  at  the  same  time  are  consistent  with  the 
sphere  solution.  Finally,  we  present  numerical  results  that  utilize  our 
pseudosymmetrlc  expansion  and  Illustrate  the  capability  of  determining  SEM 
quantities  by  patch  zoning  the  MFIE. 


i 


CONTENTS 


I.  Introduction 


II.  Pseudosymmetric  Eigenmode  solution  to  the  Magnetic 
Field  Integral  Equation 


III.  Consistency  of  Zero  Residue  at  Interior  Resonances 
and  the  Pseudosymmetric  Eigenmode  Solution 


IV.  MFIE  Pseudosymmetric  Eigenmode  Solution  to  the 
Sphere 


V.  SEM  Coupling  Coefficients  for  the  Sphere  Via  the 
Pseudosymmetric  Eigenmode  Solution  to  the  MFIE 


VI.  Eigenmode  Solution  to  the  Electric  Field 
Interodif ferential  Equation 


A Set  of  Assumptions  for  the  Generalization  of 
the  SEM  Sphere  Solution  to  an  Arbitrarily  Shaped 
Conducting  Body 


VIII.  Numerical  Determination  of  Pole  Locations  and 
Eigenfunctions  for  the  Sphere 


Appendix  A.  Proof  of  Equation  (14) 


Aooession  For 


NTIS  GPjA&I 
DOC  TAB 
Unannounced 
Justification 


Ayr  1 1 f.fclUty  Code 

Avail  and/or 

l at  special 


I 


SECTION  I 

INTRODUCTION  AND  SUMMARY 

The  original  objective  of  this  study  was  to  determine 
whether  quantities  required  by  suggested  singularity  expansion 
method  (SEM)  expansions  could  be  obtained  through  numerical  pro- 
cedures that  include  patch  zoning  the  magnetic  field  integral 
equation  (MFIE) . We  found  that  this  could  be  accomplished  at 
the  expense  of  very  lengthy  computer  runs.  Furthermore,  we 
were  motivated  to  study  the  underlying  theory  of  SEM  by  the  fact 
that  there  was  more  than  one  suggested  SEM  expansion. 

The  predominant  effort  to  provide  a theoretical  basis 
for  SEM  has  been  directed  toward  the  problem  of  determining 
the  current  density  induced  on  a perfectly  conducting  closed 
surface  by  an  incident  electromagnetic  wave  (i.e.,  the  EMP 
external  interaction  problem) . This  problem  is  treated  by 
studying  equations  that  result  from  taking  the  Laplace 
transform  of  Maxwell's  equations.  This  is  the  problem  area 
that  we  also  address  in  our  study  of  SEM  theory.  The  approach 
we  took  was  to  study  the  relationship  between  eigenmode 
expansion  method  (EEM)  solutions  and  SEM  expansions. 

We  begin  our  effort  by  studying  the  EEM  as  applied  to  the 
MFIE.  The  reason  for  this  is  that  EEM  and  SEM  can  be 
related  and  EEM  does  not  suffer  from  the  same  uncertainty  as 
does  SEM.  This  is  the  case  since  EEM  has  been  studied  as  a 
special  case  of  the  spectral  theory  of  operators.  This  does 
not  imply  that  EEM  should  replace  SEM,  because  the  suggested 
SEM  expansions  have  certain  advantages  over  EEM  expansion. 


The  ability  to  relate  EEM  to  SEM  was  facilitated  by 
recognizing  that  it  was  possible  to  use  special  properties 
of  the  MFIE  in  order  to  operationally  simplify  and  interpret 
the  standard  EEM  expansion.  Specifically,  the  standard  EEM 
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expansion  requires  that  eigenmodes  of  both  the  MFIE  operator 
as  well  as  its  adjoint  be  determined  since  the  MFIE  operator 
is  not  self  adjoint.  The  special  properties  of  the  MFIE 
operator  allow  exactly  the  same  EEM  expansion  to  be  obtained 
as  the  described  standard  solution  by  utilizing  only  the 
eigenmodes  of  the  original  MFIE  operator.  Because  the  MFIE 
operator  is  not  self  adjoint  and  we  are  still  able  to  obtain 
an  EEM  expansion  that  requires  only  the  eigenmodes  of  the 
original  MFIE  operator,  we  term  the  expansion  obtained,  a 
pseudosymmetric  eigenmode  expansion  for  the  MFIE.  The  pro- 
cedure for  grouping  and  then  utilizing  these  eigenmodes  are, 
however,  considerably  different  from  the  procedure  that  would 
be  employed  for  a symmetric  operator.  The  standard  EEM 
solution  and  the  pseudosymmetric  EEM  solution  are  always 
termwise  identical  if  explicitly  evaluated.  In  fact,  the 
eigenmodes  used  in  both  EEM  expansions  are  theoretically 
identical,  independent  of  evaluation.  It  is  the  expansion 
coefficients  that  are  operationally  different  (i.e.,  different 
representations  of  the  same  quantities) . It  is  this  difference 
that  allows  the  pseudosymmetric  representation  to  more  readily 
yield  a result  that  is  necessary  to  relate  EEM  to  SEM. 

In  order  to  discuss  this  result,  we  will  describe  the 
three  important  classes  of  quantities  that  are  required  by 
any  of  the  suggested  SEM  expansions.  They  are  the  natural 
frequencies,  the  natural  modes,  and  the  coupling  coefficients. 
The  natural  frequencies  are  the  values  of  the  complex 
Laplace  frequency,  y,  for  which  the  eigenvalues  of  the  MFIE 
operator,  A^,  is  zero.  The  natural  modes  are  the  EEM 
eigenmodes  evaluated  at  the  natural  frequencies.  The 
question  of  just  what  is  a coupling  coefficient  is  an  open 
question  that  will  shortly  be  addressed  in  more  detail.  The 
necessary  SEM  result  is  related  to  the  fact  that  Ai  has  ^eros 
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corresponding  to  both  the  interior,  Y^j>  and  the  exterior, 

closed  surface  problem  and  always  appears  in  the  denom- 
inator of  either  form  of  the  EEM  expansion  coefficients.  It  is 
essential  for  SEM  purposes  that  the  potential  singularity  in 

the  EEM  expansion  coefficient  at  yj\  be  eliminated  from  the 
final  expansion.  The  pseudosymmetric  representation  of  the 
expansion  coefficient  is  such  that  the  numerator  in  that 
representation  is  readily  shown  to  vanish  at  the  yTj  thus 
permitting  a cancellation  of  the  unwanted  singularity  and 
this  is  the  described  necessary  result.  Marin  and  Latham 
(ref.  1)  in  their  SEM  investigation  also  addressed  the  interior 
resonance  issue  and  obtained  the  same  result. 


We  now  address  the  coupling  coefficient  issue  according 
to  the  principle  that  EEM  expansions  and  SEM  expansions  should 
be  identical  in  the  Laplace  complex  frequency  domain.  The 
fact  that  we  are  eventually  interested  in  the  time  domain 
response  corresponding  to  the  inverse  Laplace  transforms, 
provides  no  theoretical  justification  for  two  different 
complex  frequency  domain  solutions  of  the  sime  induced 
surface  current  density.  There  is  essentially  a one  to  one 
correspondence  between  transform  pairs  (i.e.,  the  real 
function  of  time  and  its  Laplace  transform) . The  inverse 
transform  of  a given  Laplace  transform  function  can  result 
in  time  functions  that  differ  only  on  a set  of  measure  zero 
in  the  time  domain  (e.g.,  to  correspond  to  the  uncertainty 
of  the  function  value  at  a jump  discontinuity).  Two 
different  Laplace  transform  functions  must  necessarily  result 
in  two  different  inverse  Laplace  transform  time  signals. 


1.  Marin,  L.  and  R.  W.  Latham,  Analytical  Properties  of  the 
Field  Scattered  by  a Perfectly  Conducting  Finite  Body, 
Interaction  Note  92,  Air  Force  Weapons  Laboratory,  1972. 
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This  leads  to  two  questions  that  must  be  addressed  concerning 
the  relationship  between  the  EEM  expansion  and  suggested  SEM 
expansions.  Are  they  identical?  If  not,  is  the  difference 
significant? 

This  effort  is  primarily  directed  toward  the  first 
question.  We  address  this  question  by  obtaining  explicit 
EEM  results  and  SEM  results  for  the  case  where  the  closed 
surface  is  a sphere.  These  explicit  results  for  EEM  are 
obtained  by  using  the  pseudosymmetric  expansion  coefficient. 

The  SEM  representations  explicitly  evaluated  are  those 
employing  class  1 and  class  2 coupling  coefficients.  Obtain- 
ing the  EEM  results  was  facilitated  by  employing  the  explicit 
eigenmodes  and  eigenvalues  for  the  sphere  problem  presented 
by  Marin  (ref.  2) . The  benefits  of  the  pseudosymmetric 
expansion  coefficient  representation  became  apparent  in  this 
sphere  calculation  as  it  allowed  us  to  use  only  a limited 
portion  of  the  detailed  sphere  results  obtained  by  Marin 
(ref.  2).  The  resulting  pseudosymmetric  EEM  expansion  was 
in  exact  agreement  with  known  sphere  results.  Specifically, 
it  duplicated  the  standard  Mie  solution  and  Baum's  original 
SEM  results  (ref.  3).  This  duplication  of  results  can  only  be 
seen  after  each  representation  is  rearranged,  but  it  is  only 
rearrangement  that  should  be  permitted  in  the  theoretical  com- 
parison of  solutions.  The  agreement  of  the  pseudosymmetric 
EEM  solution  accomplished  two  purposes.  It  promoted  confidence 
in  the  pseudosymmetric  EEM  theory.  Secondly,  it  enabled  a 


2. 


3. 


Marin,  L. , Natural-Mode,  Representation  of  Transient 
Scattering  from  Rotationally  Symmetric,  Perfectly  Conduct- 
ing Bodies  and  Numerical  Results  for  a Prolate  Spheroid, 
Interaction  Note  119,  Air  Force  Weapons  Laboratory,  1972, 
Baum,  C.  E. , On  the  Singularity  Expansion  Method  for  the 
Solution  of  Electromagnetic  Interaction  Problems,  Interaction 
Note  88,  Air  Force  Weapons  Laboratory,  1971. 
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direct  evaluation  of  SEM  class  1 and  class  2 coupling  co- 
efficients (according  to  their  pseudosymmetric  description) . 

We  found  class  1 to  yield  an  exact  rearrangement  of  the  known 
solution  while  class  2 led  to  an  erroneous  expansion.  To 
conclude  that  the  class  expansion  is  erroneous,  we  argue  that 
each  sum,  corresponding  to  the  SEM  expansion  employing  either 
class  1 or  class  2 coupling  coefficients,  must  necessarily 
yield  different  results  in  the  limit  as  more  terms  are  added. 

This  can  be  concluded  without  even  examining  the  limiting 
process  since  each  class  of  coupling  coefficients  leads  to 
a termwise  different  expansion  coefficient  multiplying  the 
exact  same  linearly  independent  function.  (See  end  of  section  V.) 

Having  found  that  class  2 coupling  coefficients  based 
on  the  MFIE  gave  theoretically  incorrect  results,  we  decided 
to  examine  coupling  coefficients  for  the  sphere  problem  based 
on  the  electric  field  integro-dif ferential  equation  (EFIDE) . 

Again  we  found  class  1 to  yield  an  exact  rearrangement  of  the 
known  solution  while  class  2 led  to  an  erroneous  expansion. 

Our  investigation  of  the  sphere  problem  allowed  us  to  arrive 
at  a set  of  SEM  assumptions,  to  be  further  investigated, 
that  are  applicable  to  a general  class  of  closed  surfaces 
and  are  not  in  conflict  with  the  known  sphere  solution. 

At  this  point  we  can  conclude  that  SEM  expansions  based 
on  class  2 coupling  coefficients  are  not  identical  with  EEM 
expansions  while  SEM  expansions  employing  class  1 coupling 
coefficients  are  not  in  conflict  with  the  general  assumptions 
that  we  have  identified.  The  question  of  whether  the  theo- 
retical error  caused  by  employing  class  2 coupling  coefficients 
is  significant  was  not  specifically  addressed  in  this  report. 

In  relation  to  this  issue,  it  is  known  that  the  basic  SEM 
representation  is  of  questionable  value  for  high  frequencies/ 
early  times.  This  early  time  concern  is  further  enhanced 
because  the  time  dependent  behavior  of  each  term,  resulting 
from  the  inverse  Laplace  transform  of  an  SEM  expansion 
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employing  class  1 coupling  coefficients,  would  exhibit  a 
propagating  turn  on  time  faster  than  the  speed  of  light. 

Whether  this  deficiency  is  important  is  related  to  whether 
for  these  early  times  the  total  SEM  expansion  is  of  utility. 

A possible  rationale  for  using  the  theoretically  wrong 
class  2 coupling  coefficients  is  that  they  do  not  result  in 
this  termwise  propagation  speed  problem.  As  a consequence 
the  corresponding  SEM  expansions  have  the  potential  to  yield 
better  approximations  for  early  time  applications  than  theo- 
retically correct  expansions.  The  works  of  Tesche  (ref.  4) 
and  Marin  (ref.  2)  were  not  specifically  directed  at  this 
issue;  however,  they  provide  some  evidence  that  SEM  expansions 
utilizing  class  2 coupling  coefficients  potentially  can  yield 
meaningful  time  domain  solutions.  Both  references  considered 
time  domain  solutions  for  the  total  current  induced  on  an 
object,  a thin  wire  in  (ref.  4)  and  a prolate  spheroid  in 
(ref.  2).  In  either  case  they  compared  the  time  domain 
solutions  obtained  by  an  alternate  procedure  to  the  solution 
obtained  by  an  SEM  expansion  employing  class  2 coupling 
coefficients  and  they  obtained  good  agreement.  At  this  point, 
the  use  of  class  2 coupling  coefficients  should  be  viewed 
as  an  empirical  approximation. 

We  view  the  status  of  SEM  as  providing  suggested  repre- 
sentations having  potential  benefits  to  a variety  of  wave  inter- 
action applications  as  well  as  to  electromagnetic  pulse  (EMP) 
interaction  problems;  however,  there  is  at  present,  insufficient 
knowledge  to  decide  a priori  how  well  any  suggested  SEM  expan- 
sion represents  even  external  interaction  quantities. 


4.  Tesche,  F.  M. , On  the  Singularity  Expansion  Method  as 
Applied  to  Electromagnetic  Scattering  from  Thin  Wires, 
Interaction  Note  102,  Air  Force  Weapons  Laboratory,  1972. 
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SECTION  II 


PSEUDOSYMMETRIC  EIGENMODE  SOLOMON  TO  THE  MAGNETIC 
FIELD  INTEGRAL  EQUATION 

In  this  section  we  derive  an  eigenfunction  expansion  solution 
to  the  Magnetic  Field  Integral  Equation  (MFIE)  by  employing  a 
procedure  that  allows  the  determination  of  the  expansion 
coefficients  without  explicitly  calculating  the  eigenfunctions 
of  the  adjoint  operator.  We  call  this  solution  the  pseudo- 
symmetric  eigenmode  expansion  for  the  MFIE. 


We  start  with  the  MFIE  for  the  exterior  problem 


j J(r)  - i n (r)  * [vG(r,r')  * J (r ' )]  dS • = Jlnc(r)  (1) 

•'s 

where  J(r)  is  the  induced  surface  current  density  on  the  perfectly 
conducting  surface  S,  n(r)  is  the  outward  unit  normal  of  S, 

Jinc  = n * Hlnc,  Hlnc  is  the  incident  magnetic  field  and  G is 
the  free  space  Green's  function, 


G(r,r') 


1 e-Y|r-r' 

4 7T  | r - r 1 | 


We  rewrite  equation  (1)  as 


5?J  = 


jinc 


(2) 


with 

<f  = \ - L (3) 

where  I is  the  identity  operator.  The  solution  to  equation  (2)  is 
given  by 


L— j 


where 


TJ.  = A.J.  (5) 

—i  l—i  ' ' 

The  expansion  coefficients  a.  involve  an  inner  product  for  Jlnc 

1 + 

and  the  eigenfunctions  of  the  adjoint  operator  !f  . As  we 
mentioned  in  the  beginning  of  the  section  we  need  not  determine 
these  eigenfunctions  explicitly  because  we  can  employ  the  pseudo- 
symmetric  method  which  we  will  now  outline.  We  begin  by  defining 
operators  M and  Q through 


Mf  = n x f 

Q = ML 


(6) 


2 /v  ~ 

where  f_  is  a surface  vector.  Noting  that  M£=nxnxf^=-f^ 
leads  to  the  following  property 


M 


2 


-I 


(7) 


Next  we  operate  on  equation  (5)  with  M and  obtain  the  eigenvalue 
equation 


M = XiMJi 
M.y.'=  | M - 


(8) 


In  conjunction  with  equation  (8)  we  define  the  adjoint  equation 


(M*')+F. 


(9) 


where  (MV')  , M:  obey  the  usual  adjointness  relationship 


(f,Aa)  = (ATf,2) 


(10) 


with  an  inner  product  defined  as 


(a ,b)  = J a* 


• b dS 


With  the  aid  of  equations  (8),  (9),  and  (10)  we  can  derive  the 
relationship 


(*i  - Xj  ) (Lj/MJ^)  = 0 


(11) 


For  generality  we  assume  that  the  eigenvalues  are  degenerate  and 
rewrite  equation  (11)  as 


•!** 


ui  - Xj  > - 0 


where  Z,  m signify  degeneracy.  Using  this  one  can  show  that 


(12) 


Ai  “ Ai 


(F . „ ,MJ . ) = N . „ 6 . . 6. 
—jZ  — lm  jZ  ji  fm 


(13) 


where  the  biorthogonality  relationship  for  i = j,  Z f m is 
obtained  via  a Gram-Schmidt  biorthogonalization  procedure.  We  are 
now  in  a position  to  derive  an  important  result  by  noting  that 


(M 


'/')  1 = (|  M - o) ' = | M ' " Q ' = 


- ± M - Q 


(14) 


(see  appendix  A)  and  invoking  the  first  of  equations  (13)  to  rewrite 
equation  (9)  as 
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(15) 


(- 


Ei 


l 


MF. 

— l 


l 


By  first  operating  on  equation  (15)  with  M and  then  complex 
conjugating  the  resulting  equation,  we  obtain 


(I + - xi  eI, 


(16) 


where  we  have  used  equations  (6)  and  (7) . Equation  (16)  can  be 
rewritten  as 


(1 


- Xi^i* 


(17) 


(18) 


where  = 1/2  - L is  the  MFIE  operator  for  the  exterior  problem. 

^ ★ , 

Thus  J . 0 2 F.„  is  an  eigenfunction  of  the  MFIE  operator  with  an 

1 36  1 36 

eigenvalue  1 - where  X^  corresponds  to  J^m.  According  to  the 
second  of  equations  (13) 

* Jl ‘ * iim  ds 

s 

- fln  ■ S * iim  dS  - NU  5lj  S£m  (19> 

and  we  define  this  integral  form  as  a pseudo  inner  product 


n * j . dS 
— lm 


N 


il 


(20) 


We  now  return  to  equation  (2)  to  obtain  the  eigenfunction  expansion 
solution . 


as  the  pseudosynunetric  eigenmode  expansion  solution  to  the  MFIE 
or  equivalently  to  equation  (2).  The  eigenmodes  required  by  the 
expansion  are  eigenmodes  of  the  original  operator  W taken  in 


pairs,  one  corresponding  to  the  eigenvalue  and  the  other  cor 
responding  to  the  eigenvalue  1-A . . 


>A 


SECTION  III 

CONSISTENCY  OF  ZERO  RESIDUE  AT  INTERIOR  RESONANCES  AND  THE 
PSEUDOSYMMETRIC  EIGENMODE  SOLUTION 


1 


r 

y 


y \ 

i 


In  this  section  we  show  that  the  set  of  eigenmodes  for  the 
exterior  problem  is  identical  to  the  set  of  eigenmodes  for  the 
interior  problem,  the  eigenvalues  A^  of  the  exterior  eigenvalue 
problem  have  zeros  at  the  inttn-Lofi  resonances  (in  addition  to  the 
exterior  resonances)  and  that  the  coefficients  in  the 

pseudosymmetric  eigenmode  solution  given  by  equation  (24)  are 
zero  at  these  interior  resonances.  The  last  result  is  necessary 
in  order  to  relate  the  EEM  to  the  SEM.  Thus  the  eigenmode  series 
may  be  rearranged  to  be  written  as  a singularity  expansion  (which 
is  a special  case  of  the  Mittag-Lef f ler  theorem  stated  in  section 
VII,  equation  (103))  which  will  not  involve  interior  resonances  since 
the  excitation  coefficients  {£^£/JinC}  will  be  shown  to  be  zero. 

We  start  by  recalling  equations  (5) , (16) , (18)  and  the 
eigenvalue  equation  for  the  interior  problem. 


"i  5 (§  - L)ii  * Xi-i 

(I  + L)ii  - xiii 
xzil  = (f  + l)jx  - 


(5) 


(25) 


(26) 


where  the  superscript  "I"  denotes  "interior."  To  simplify  the 
notation  we  will  not  supply  a superscript  "E"  to  exterior  quanti- 
ties and  also  ignore  degeneracy  because  our  results  do  not  depend 
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on  degeneracy.  The  eigenvalue  X^  in  equation  (5)  is  any  eigen- 
value of  the  exterior  operator  and  either  the  ' s or  ' s comprise 
the  entire  set  of  eigenmodes  of  the  exterior  operator  rJt . Bearing 
this  in  mind,  a comparison  of  equations  (25)  and  (26)  shows  that 
the  interior  and  exterior  eigenvalues  comprise  the  same  set  and 
that  the  interior  and  exterior  eigenfunctions  comprise  the  same 
set. 


We  now  proceed  to  show  that  the  excitation  coefficient  evalu- 
ated at  the  interior  resonances  is  zero.  The  interior  resonances 
X?j  are  solutions  to  the  equation 

X.ty1.)  = 0 
l 1 lj 

and  the  corresponding  interior  natural  modes  satisfy  equation  (26) 
with  ^ = 0 or  equivalently  equation  (25)  with  X^  = 0. 

V"1  (Y  J j ) Ji(Y*j)  = 0 (27) 

The  solution  to  equation  (27)  can,  using  an  appropriate  normaliza- 
tion, be  given  by 


2i<Y*j,r>  - ij.  x 


-n  x Hi(yij,r)  , r e S 


(28) 


where  we  have  simply  made  the  distinction 


n 


I 


-n 


and  n is  the  outward  normal  to  S. 

H.(Y*->r),  where  r € V and  V is  the 
— i i j — — 


The  equation  satisfied  by 
interior  region  bounded 


Vxv  x HA ( y ^ j , r ) + (Y^j)2  H± ( y ^ j » r ) = 0 , r £ V (29a) 


nx  x ^ x Hi(Y^,r)j  = 0 , r £ S . 


(29b) 


The  excitation  coefficient  in  equation  (23)  evaluated  at  y^  is 


JinC  (Y 


Ji"  /-i(YijJ  ' " * iin°lYjj 


) dS 


*/S'[~ i'Yij’  * SlnC<vf  j>]  dS 


* 


(30a) 


(30b) 


In  order  to  show  that  {y^j)  — 0 we  employ  the  following  identity 
7 *[h.  x (7  x ElnC)  - ElnC  x (V  x H± )J 

= M-i  jv  * V x EinC  + (Y^j)2  ElnCJ 

- ElnC  • [v  x v x h.  + (Y^j)2  hJ  rev  (31) 

1 nC 

Since  the  source  for  the  incident  electric  field  E lies  out- 
side V and  satisfies  equation  (29a) , it  follows  that  an 
integral  of  the  righthand  side  of  equation  (31)  over  V yields 
zero.  The  volume  integral  of  the  lefthand  side  of  equation  (31) 
can  be  converted  to  a surface  integral  using  the  divergence 
theorem  to  obtain 

fn  . [hjl  x (v  x EinC)  - ElnC  x (v  x H i ) ] dS  = 0 
S 1 

inc  inc 

Substituting  the  Maxwell  equation  V x e = - YZQH  into  this 
relationship,  we  obtain 


(Y-->  = 0 
1 ij 


into  this 


SECTION  IV 


MFIE  PSEUDOSYMMETRIC  EIGENMODE  SOLUTION  TO  THE  SPHERE 


In  this  section  we  apply  the  pseudosymmetric  eigenmode  solution 
developed  in  Section  II  to  obtain  the  induced  surface  current  on 
the  surface  of  a perfectly  conducting  sphere  illuminated  by  a 
plane  wave  that  corresponds  to  the  Laplace  transform  of  a delta 
function  plane  wave  pulse.  The  pseudosymmetric  approach  produces 
the  Mie  solution  in  the  form  given  by  equation  (B-38)  in  reference  3. 


We  repeat  here  some  of  the  equations  in  Section  II  needed 
in  this  section. 

The  solution  to  the  MFIE 


where 


V'J  = J 


me 


(2) 


is 


i?  -rinci 

* ■ EE-Srvr1  fe 

iS,  1 


(23) 


*iil  “ Ai^i* 


*Zu  = (1  - 


(5) 

(17) 


5 /?«  ' " x -jm  ds 


= N. „ <S . . 6„ 
iS,  jlj  £m 


(18) 


/~  * _inc 

-it  * n - 


dS 


(33) 
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i 

> 


* 

l 

V 


I 


An  interesting  property  for  the  sphere  that  will  be  used  shortly 


where  er  s n.  Next  we  observe  that  because  of  the  Wronskian 
relationship 


u„u>[unu>]'  - un(t)  [an(?)]’  - 


the  following  is  true 


XR  + XQ  = 1 
n n 


Thus  according  to  equations  (17a)  and  (38) 

'£ R = (l  - XR)  R = XQ  R (3S 

— n,m,o  \ n/  — n,m,c  n — n,m,a 

<£Q  = (l  - XQ)  Q = XR  Q (4( 

^n , n f a \ n / — n , n , o n — n , n , a 

Comparing  equations  (39)  and  (40)  to  equations  (36)  we  conclude 


that 


* * 

— n,m,a  **  cm',a'  — n,m',a' 


m’=0  a 1 =o ,e 


— n,m,o  ^ dm',a'  — n,m',a' 


m ' =0  a ' =o , e 


We  next  show  that  the  only  nonzero  coefficients  in  equation  (41) 

are  c and  d . This  can  be  accomplished  by  first  recalling 
in  f o in  f o 

equations  (18) 


j ^n,m,o  * ^er  X -n,m",o"^  dS 


N 6 „ 6 

n,m,o  m,m"  a ,o' 


/O  • (e  x Q ..  .. ) dS 

^n(m,o  r — n,m",on/ 


Nu  6 „ 6 „ 

n,m,o  m,m"  aro" 


which  with  the  aid  of  equations  (37)  become 


* £n,m" ,a"  dS  Nn,m,o  6m,m"  6o,o" 


n'm'°  * — n,m"  ,o" 


dS  = -Nu  5 „ 5 

n,m,o  m,m"  0,0'* 


Substituting  equations  (41)  into  equations  (43)  and  employing  the 
orthogonality  properties 

y*-n, m'(a'  * -n,m",o"  dS  = * Sn,m",o"  dS 


Mn/m"o"  ^m'^m’  ^a",o' 


given  by  equations  (B-19)  in  reference  3,  we  obtain 


cm",a"  Mn,m",a"  ^n,m,o  ^m,m"  ^o,a" 


dm" , a"  Mn,m",a"  ”Nn,m,a  6a,a" 


which  prove  that  only  c and  dm  are  nonzero.  Thus  we  can  replace 

Hi  f v IR  f 0 *■ 

(41)  by 


R = dn  Q 

— n , m , a R — n , m , a 


^n,m,o  dQ  — n,m,a 


where  dR,  d^  are  arbitrary  constants. 

R 0 

Next  we  calculate  Nn,m,0,  N“  a (which  appear  in  expansion  equation 
(23))  by  invoking  equations  (45)  and  (46) 
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n,m,a  m,a  n,m,a  R n,m,a 


(47! 


N 


n,mf  a 


-d  M = -d^  M 

mfa  n,m,o  Q n,m,a 


In  order  to  calculate  the  expansion  coefficients  in  equation  (23) 
we  first  recall  equation  (33)  and  employ  equation  (46) 


{ R m JlnC> 
n / m f o ■ 


= [r 
J — n,m ,a 

- -dR  fa.m.o  ■ a1 


(er  x jlnc)  dS 


inc 


dS 


since 


* inc  A ^ inc 

e_  x = e x e x h 

r — r r — 


-Hlnc  + er(er  • Hlnc) 


and  Q is  a surface  vector.  Similarly, 

— n,m,o  J 


2n,m,o'  ilnC>  = -dQ  fZn.m.o  ' ^ 


{Q 


inc 


dS 


(48! 


(49 


The  form  of  the  incident  magnetic  field  can  be  found  by  using 
equations  (B-55)  and  (B-58)  in  reference  3 


co  n 

■EE  £ [■ 


l2  ,n ' ,m' 


M(}> 
- P -n ' , 


m 


n=l  m'=0 


- A N(1)  1 

l,n'  ,m'  ,o'  ,p  -n'  ,m'  ,o'J 


(50 


where  p = 2,3  represents  the  two  possible  polarization  directions 
defined  in  reference  3 and 


I * 
I > 


< 

L , ? 

K J 

l ' •’ 
\ ' : /\> 

' 


M(1)  „ = i (yr)  R (0,<J>) 

— n / m / o n ' — n , m , o 


(1)  in(Yr) 

; = n (n  + 1)  — P (0,d>) 

— n,m,o  yr  — n,m,o 


-n(Yr>] 


— Q ( 0 , d> ) 

yr  —11,111,0 


-n  m o is  the  third  vector  spherical  harmonic  (the  other  two  being 
Q and  R;  see  equations  (R-ll) , (B-12) , and  (B-13)  in  reference  3). 

Employing  the  orthogonality  properties  among  P,  Q,  and  R 
given  by  equations  (B-18  and  (B-19)  in  reference  3 and  equations  (50) 
and  (51),  we  can  rewrite  equations  (48)  and  (49)  as 


[Yaln(yal1 


^ £>  = d A M n 

— n,m,o'  -p  R l,n,m,o,p  n,m,o  ya 


r ~ inc  i 

^n,m,o'  — p ~ dQ  A2,n,m,o,p  Mn,m,o  in^Ya^ 


where  a is  the  radius  of  the  sphere.  We  are  now  in  a position  to 
rewrite  equation  (23)  in  its  final  form  by  first  recalling 
equations  (47)  and  (52) 


n=l  m=0  o=e , 0 


A.  _ d_  M |~yai  (ya)] 

l,n,m,o,p  R n,m,oL  n 1 J 


ya  A d M 

n R n,m,o 


R 

— n,m, o 


A2,n,m,o,p  dQ  Mn,m,o  an^Ya^ 

A°(-d^  M ) 

n Q n,m,o' 


and  using  equation  (35)  to  obtain 


Q 

— n , m , o 


Jp(e,<t>) 


EE  E 

r»  = l m=  n rr  — ea 


n, 111,0, p 
2 kn(ya) 


R m „(0,<f>) 
— n t in  t o 


1 

2,n,m,o ,p  0 ( 

i / \V  — n,m,o  (0  ,<t>)  < 
ya  ya  (ya)  ' ,y  ) 


(53) 


Equation  (53)  is  identical  to  equation  (B-68)  in  reference  3. 
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SECTION  V 

SEM  COUPLING  COEFFICIENTS  FOR  THE  SPHERE  VIA  THE  PSEUDOSYMMETRIC 

EIGENMODE  SOLUTION  TO  THE  MFIE 

We  first  examine  Class  1 coupling  coefficients  with  q = 1 

corresponding  to  the  TM  modes  R . In  appendix  B we  show  that 

— n , m , o 

this  coupling  coefficient  has  the  form 


(v) 

n,n ' ,m,o ,p  ' 


(Y^.n,-Y)ctn( ^~n,m,o’  -p 

|nR  [dAR/dyl  ) 

1 n»m,a  [_  n/  TJ  ) R 

Y Yn,n' 


where  t = -a/c.  From  the  previous  section  we  recall  equations 
(52),  (47)  and  (35) 


(R  . J } = A..  d_  M — 

— n,m,o  — p l,n,m,a,p  R n,m,a  ya 


N m „ = d_  m _ (independent  of  y) 
n,ra,o  R n,m,a  r 

An  = (W^ya)]'  [ya  kn<ya>] 

Recalling  that  for  q = 1 modes  ya  k (ya)  R = 0 we  have 

L n JY  Y.  ni 


[dAn/dY]y=yR>n,  = {a[Yain(Ya)J’  *n(Ya)] 

+ a|^yain(ya)j  j^ya  kn(ya)J  jY=yR 

■ a{[Yain<H’  [va  kn(ra!]jy-y®,n. 
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In  view  of  equations  (55),  (56),  and  (58),  equation  (54)  gives 


RU>  M = e'Vn^’^o 

nn,n'  ,m,o ,plY 


A /a 

l,n,m,o ,p 


|ya[ya  kn(yal] 


R 

Yn,n’ 


Following  equation  (B-76)  in  reference  3 in  conjunction  with 
equation  (B-74)  we  find 


(ya)2  kn(ya)  = e"Ya  C1>n(ya) 


Noting  that 


(ya)  kn(ya) 


fn,n' 


= < ya  kn(ya)  + ya|ya  kn(ya)j 


'n , n ' 


= jyajya  kn(ya)j 


'n , n ' 


we  can  rewrite  equation  (60)  as 


1 i ,1 

(ya[ya  VYa)]’  LyR 

' Tn,n' 


'(Ya)2  kn(ya)]’  [e'Ya  Clin(ya)]' 

• • y*  r-\  • 


nf  n 1 


K,a[cl,n(Ya>J' 


- e‘Y“  C. 


„(Ya)) 

,n  j 


R 

Y=Yn, n’ 


Combining  equations  (59)  and  (61)  we  obtain 


„ n _(Y)  = e 
n,n  f rn / o t p 


(Y*  _ i -Y) ct 


A.  /a 

1 , n , m , a , p' 


f „ _ | 

n,n 


A,  /a 

l,n,m,a,p 

[Cl,nlYa)]'  R 
L JY=Yn,n’ 


-yet 
e o 


since  t = -a/c, 


Employing  equation  (B-82)  in  reference  3 in  conjunction  with 
equation  (B-98)  we  finally  obtain 


R (1)  . . _ -yet  0 

nn, n ' ,m, o ,p  Y e ° 


l,n,n'  ,m,a,p  c 


Recalling  that  y = s/c  we  understand  that  the  coupling  coefficient 
given  by  equation  (63)  is  identical  to  the  one  that  can  be  identi- 
fied in  equation  (B-97)  in  reference  3 for  q = 1. 

Following  the  same  procedure  as  above  we  obtain 


n°  , (Y)  - e‘,cto  c,  , i 

n,n',m,a,p  2,n,n’,m,o/p  c 


Thus  Class  1 coefficients  give  the  correct  result  for  the  sphere, 
i.e.,  equation  (B-97)  in  reference  3.  For  Class  2 the  defining 
relationship  for  q = 1 is 


_R ( 2 ) ,yy  _ — n,m,o  * ^p  J 

nn,n'  ,m,a,plY  TZ R 
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Taking  into  account  the  manipulations  that  led  to  equation  (63) 
we  can  rewrite  equation  (65)  as 


n?«>.  . „ (Y)  = ± aY"">,a  c 

n,n  ,m,o ,p  c 


l,n,n' ,m,o ,p 


hrai^Ya)J/Ya 

{[yain(yal]'/va\ 

1 ’n,n' 


and  similarly  for  q = 2 


nQ(2)  (y) 

n,n' ,m,o,p  1 


I a n,n’ac 


in(Ya) 


2,n,n',m,o,p  Q 

n n , n 


Comparing  these  expressions  with  (63)  and  (64)  we  conclude  that 

Class  2 does  not  give  the  correct  answer  for  the  sphere.  (See  end 

of  this  section.)  Specifically  let  us  consider  q = 2,  n = 1, 

n'  = 1 and  set  Y?  = Y,  . We  have 

n / n 


ii(va)  - £n?ai  - - & f1  - &) + %£  f1  + h) 


To  Laplace  invert  the  form  containing  the  Class  1 coupling 
coefficient  for  n = 1,  n'  = 1,  we  write 


,(D 


I c eYa  -)  = C e*1 

c 2 v - Yl)  2 


c (t+a/c) 


u(t  + a/c) 


where  C-  = C„  , „ „ with  n = 1,  n’  = 1.  We  also  invert  the 

2 2,n,n  ,m , a , p 

term  n = 1,  n’  = 1 with  the  Class  2 coupling  coefficient  by 
considering  the  expression 


inUia)  1 2ya  \ yal  Y ~ Y^  2Ya  \ ya)  Y **  Yx 

(70) 


Now  we  can  write 


/a  _ _L\ 

1 

1 

(ayl 

- 1 

aYl  “ 

1 

+ Ii\ 

W Y2) 

Y 

Y1 

” 2 
Y1 

l Y - 

Y1 

Y 

y2> 

/a  + -A-\ 

1 

1 

/aYl 

+ 1 

ay1  + 

1 

- Ii) 

U Y2j 

Y~ 

*1 

' Y - 

Y1 

Y 

Y2' 

and  equation  (70)  can  be  rewritten  as 


C,  eYla 
t(2)  - ^ 2 2 


2Y1a  in(Y1a) 


YjC (t+a/c) 

(aY-j^  - l)e  + 1 + Y 


lct]u 


(t  + a/c) 


[(aYi 


+ 1)  e 


Yj^c  (t-a/c) 


- 1 - YictJu(t  - a/c) | 


fae 

2Y 


Y^c (t+a/c) 


>e  r Yla 

Y~ = aY-j^  - l)e  u(t  + a/c)  + (ay^^  + 

ia  W 


1) 


-Y!a 


1 u(t  - a/c)  + (1  + Ylct)(u(t  + a/c)  - u(t  - a/c)jj 


(71) 


For  -a/c  < t < a/c  equation  (71)  is  different  from  equation  (69) . 
However,  for  t > a/c  equation  (71)  can  be  reduced  to 


,?v  y-ic  (t+a/c) 

Tl2)  = C2  e 
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and  Class  2 gives  the  correct  answer  only  for  t > a/c,  i.e., 
after  the  wavefront  has  passed  the  sphere.  Notice  that  the  previous 
result  is  valid  for  a delta  function  incident  plane  wave.  For  an 
incident  wave  with  a different  functional  dependence  say 


fp(t  - r/c)  u(t  - e-^-r/c)  the  response  is  split  into  the  object 


response  and  waveform  response,  i.e.. 


i<>p  - £yy  vs)  ve>(s  - v 


£VS) 


f (s)  - f (s  ) 
p P a' 


s - s 


From  equation  (72)  we  see  that  Class  2 will  again  give  the  correct 
object  response  only  for  t > a/c.  However,  the  waveform  response 
will  be  wrong  for  all  times  since  from  equation  (73)  we  have  to 
convolve  & 1 (na(s)/(s  - sa))  with  1 (f  (s)  - fp(sa))  and 


(na  (s)/Ts  - sa))  is  correct  for  t > a/c  only. 


We  conclude  this  section  by  answering  the  following  question 
What  if,  despite  the  fact  that  Class  1 and  Class  2 coupling  co- 
efficients are  different,  the  corresponding  infinite  sums  produce 
identical  responses?  The  answer  is  that  this  is  impossible.  To 
show  this,  we  recall  the  orthogonality  properties  of  the  v^'  s 


nmo '-nmo^ ' i*e"  equations  (B-19)  in  reference  3 (which,  inciden- 


tally, show  the  linear  independence  of  the  v^’  s as  we  mentioned  in 
the  introduction)  and  notice  that  if  the  responses  were  identical 


we  would  have 


An  , 

Sn  / n / in  y o y | 

y-y_, 


where  An  . is  the  difference  between  the  Class  1 and 

'n,n'm,o,p 

Class  2 coupling  coefficients.  In  particular,  for  q = 1 

(R  functions)  and  n = 1,  we  have  n'  = 0;  i.e.,  there  is  only 
nmo 

one  pole  and  it  lies  on  the  real  axis.  Thus  An?  n m „ ~ 0. 
c 1,0, m,o ,p 

This  is  impossible  as  we  have  shown,  and  consequently  the  two 
sums  must  be  different. 


SECTION  VI 
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EIGENMODE  SOLUTION  TO  THE  ELECTRIC  FIELD 
INTERODIFFERENTIAL  EQUATION 

The  Electric  Field  Integrodif ferential  Equation  (EFIDE)  for 
scattering  from  perfectly  conducting  bodies  can  be  cast  into  the 
following  operator  form  (ref.  2) 

Z • J = z"1  E*nc  (74) 


where 


Z • J » I • (V4>  + yA) 


A 


I - nn 

j G (r , r ' ) J (r 1 ) dS' 
S 


4> 


7 / G(r,r') 
S 


V’ 


J(r')  dS’ 


inc  inc  • • 

and  Et  is  the  tangential  (=  E ) component  of  the  incident 

electric  field.  It  can  be  shown  that  the  following  relationship 

is  true 


(75) 


Thus  if  we  define  a complex  inner  product 


we  have 


(a, Z • b)  = (Z*  • a,b) 
32 


(76) 


* 


( 


•- 

-« 


In  order  to  utilize  equation  (76)  in  the  derivation  of  the  eigenmode 
solution  to  equation  (74)  we  define  the  eigenvalue  equations 


S ' in  * ci  in 


2+  • 4*  * 4 in 


and  one  can  show  that 


t* 

«i  = «i 


( t* 

i * ^jm  dS  " NU  6ij 


Recalling  equation  (76)  and  the  defining  relationship 


(a,Z  • b)  = (Z  • a,b) 


we  understand  that 


Z+  = Z* 


We  can  now  rewrite  equation  (78)  as 


Thus  in  view  of  equations  (79)  and  (80)  the  eigenmode  solution 
to  equation  (74)  is 


where 


(8: 


(8: 


As  an  example  for  equation  (81)  we  consider  the  case  of  a sphere 
From  reference  2 we  have 

n = 1,2,  ...»  00  I 

m _<  n / (8 

o = o,e  (odd,  even) I 


j ^n,m,< 

1° 

\ — n,m,< 


i.e.,  J!  has  the  same  eigenfunctions  as  the  MFIE  operator  fJT  . The 
corresponding  eigenvalues  are  (ref. 2) 


’ [Yain(Ya)]  [Ya  kn(Ya)] 

cn  = ‘[Yatn(Yal]  [Ya  VYa>] 


(8 


and  from  equation  (44) 


Nifc  Mn,m,a  J"-n,m,a  ’ ^n,m,< 


dS 


= I Q • Q dS 

J ,m, a — n , m , a 


(8 
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The  incident  electric  field  Einc  has  the  form  (equation  (B59)  in 


ref. 2) 


Einc=yy>  V [a.  „ M^>  (yr) 

— Z— i L-J  Z— / l l,n,m,0  ,p  — n,m,a  ‘ 


n=l  m=0  a=e , o 


+ A-  N(1)  (yr)l 

2,n,mfa,p  — n,m,o'  — J 


where  and  are  given  by  equation  (51)  and  the  A's 

by  equation  (B-58)  in  reference  3..  By  using  the  orthogonality 
relationships  given  by  equations  (B-ll) , (B-12) , and  (B-13)  in 
reference  3 in  conjunction  with  equation  (86)  we  find 


C<m,o'  ^nC>  = V^n^o'  ^ 


in(ya)  Mn>mf0  Ai/n,m,0,p 


C^n,m,a'  ^ = C<m,0'  ^ 


(87a) 


_ [^(Ta)] 


M A0 

n,m,o  2,n,m,o,p 


(87b) 


where  M are  given  by  equation  (85) . Employing  equations  (87) 

n ,m,a 

and  (84)  we  can  rewrite  equation  (81)  as 


= y'  y in<,<al  r 

^ “i  M ("vai  (ya)l  [ya  k (ya)l  ~n ' m ' 0 

n=l  m=l  o=o, el  n,m,a[_  n J [ n' 

- A2 ,n,m,a,p{[7a:i'n  ^ a } Mn,m,o  1 


’.....K'’'11]!''1  *„<**>]' 


V V V'  l , n ,m, a , p 

t-J  ,va)2  k (va)  -n,m,o 

n=l  m=l  o=o, e L(Ya}  knlY  ' 


A„ 

2 , n , m , o , p 

T , , XV  — n,m,o 

j^ya  kn(ya)  I 


Equation  (88)  as  expected,  is  identical  to  equation  (53)  obtained 
via  the  pseudosymmetric  eigenmode  expansion  for  the  Magnetic  Field 
Integral  Equation. 

We  conclude  this  section  by  noting  two  important  factors; 

(a)  if  we  examine  the  coupling  coefficients  for  the  sphere  we 

can  follow  a procedure  similar  to  the  one  employed  in  ?.  ,ction  IV 

and  arrive  at  formulas  that  are  identical  to  the  ones  in  section  IV, 

i.e.,  draw  the  same  conclusions  and  (b)  the  excitation  coefficient 
* “*1  inc 

(J i,  Zo  E^  ) can  be  shown  to  vanish  at  the  interior  resonances  by 
employing  a procedure  similar  to  the  one  in  section  III  of  the 
MFIE. 
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A SET  OF  ASSUMPTIONS  FOR  THE  GENERALIZATION  OF  TIIE  SEN  SPHERE 
SOLUTION  TO  AN  ARBITRARILY  SHAPED  CONDUCTING  BODY 

In  this  section  we  cast  the  eigenmode  solution  (to  either  the 
MFIE  or  EFIDE)  for  electromagnetic  scattering  from  an  arbitrarily 
shaped  perfectly  conducting  body  into  a form  that  represents  a 
generalization  to  the  SEM  sphere  solution.  Certain  assumptions 
are  made  along  the  way  that  are  motivated  by  the  known  sphere 
solution  and  the  procedure  leads  to  an  SEM  representation  that 
involves  Class  1 coupling  coefficients  with  no  additional  entire 
function.  If  the  same  procedure  is  applied  to  the  sphere,  no 
assumptions  are  necessary,  and  one  is  inexorably  led  to  the 
SEM  solution  with  Class  1 coupling  coefficients  with  no  entire 
function  to  be  added. 

We  begin  with  the  eigenmode  solution 

J(r,y)  = J±  (89a) 

where  b.  = {J.,  Jinc}/N.  for  the  MFIE  (eq.  23)  and  b.  = (J*,  Z_1 
inc  ^ l o 

)/N^  for  the  EFIDE  (eq.  81).  The  eigenvalues  and  eigen- 
functions are  determined  by  solving  the  appropriate  eigenvalue 
problems.  It  has  been  shown  in  reference  1 that  for  an  incident 
delta  function  plane  wave  pulse  (in  the  time  domain)  J(r,y) 
is  a meromorphic  function  of  y,  i.e.,  in  any  finite  region  of  the 
complex  yplanc  J(r,y)  has  a finite  number  of  pole  singularities. 

The  pole  locations  correspond  to  the  exterior  and  interior  resonances 
of  the  body  and  can  be  determined  by  setting  A^(y)  = 0.  Recall  that 
our  incident  wave  in  the  time  domain  is  a delta  function  plane  wave 


3 


pulse,  i.e.,  the  incident  fields  have  no  singularities  in  the 
finite  Y-plane.  In  section  III  we  showed  that  the  excitation 
coefficient  (J^,  J^nc)  is  zero  at  the  interior  resonances 
and  in  section  VI  we  mentioned  that  we  can  similarly  prove  that 


(J i»  ZD1  E"*"nc)  is  also  zero  at  y*j, 


This  leads  us  to  one  of  the 


assumptions  necessary  for  our  generalization;  it  is  that  J(r,y) 


can  be  written  as 


s • J (r,y) 


= 2 Li(- 


(89b) 


where  s is  a unit  surface  vector  and  the  L^(r,y)  are  meromorphic 
functions  of  y.  At  this  point,  it  is  appropriate  to  relate  equa- 
tion (89b)  to  the  Mittag-Lef f ler  theorem  (see  refs.  5 and  6).  This 
theorem  is  used  to  derive  a general  representation  for  meromorphic 
functions  in  terms  of  an  infinite  sum  in  which  poles  are  explicitly 
represented.  After  we  obtain  more  explicit  representations  for  the 
L.(r,y),  we  will  show  that  equation  (89b)  does  not  violate  this 

1 . T 

theorem.  Because  the  excitation  coefficients  are  zero  at  y . . , at 

i] 

this  point  we  assume  that  the  L^(r,y)  are  such  that  their  only  singu- 
larities in  the  finite  y-plane  are  poles  located  only  at  exterior 
resonances.  Notice  that  to  each  subscript  i corresponds  a set  of 

exterior  resonances  y . . . 

ID 

Equation  (89b)  is  exact  for  the  sphere  as  we  can  see  by  invoking 
equations  (34),  (35),  and  (52)  for  the  MFIE  and  equations  (83),  (84), 

and  (87)  for  the  EFIDE.  Next  we  invoke  a corollary  to  the  Weierstrass 
theorem  concerning  the  representation  of  an  entire  function  which 
states  that  "every  function  which  is  meromorphic  in  the  whole 

5.  Carrier,  G.,  M.  Krook  and  C.  Pearson,  Functions  of  a Complex 
Variable,  McGraw-Hill,  New  York,  1966. 

6.  Ahlfors,  L. , Complex  Analysis,  McGraw-Hill,  New  York,  Second 
Edition,  1966. 
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finite  plane  is  the  quotient  of  two  entire  functions"  (see  for 
example,  ref.  6).  Thus 


(r,Y) 


Li^'Y>  * g,  (VT 


(90) 


and  the  zeroes  of  g^(y)  correspond  to  the  poles  of  . Equation 
(90)  is  also  true  for  the  sphere.  We  now  assume  that  both  f^ 
and  gi  have  a finite  number  of  zeroes  for  each  i and  consequently 
being  entire  functions  they  must  have  the  form  (ref. 6) 


fi(r,y)  = eFi(^'Y)  Qi(r,Y) 


(91) 


gi(Y)  = eGi(Y)  Pi(Y) 


where  , Gi  are  entire  functions  and  polynomials.  Guided 

by  the  sphere  solution  we  assume  that  F.  is  a function  of  y only 
and  that  the  degree  of  Q is  lower  than  the  degree  of  P in  order 
to  be  able  to  obtain  a partial  fraction  expansion  of  the  desired 
form.  With  the  aid  of  equations  (91)  we  can  rewrite  equation 
(90)  as 


Li(r,Y) 


. ♦i(Y) 


Qi (r,Y) 


( Y) 


(92) 


where  we  have  defined  4k  = F^  - G^ . 


Again  equation  (92)  is  exactly 
true  for  the  sphere.  The  explicit  form  for  a sphere  can  be  obtained 
by  invoking  equations  (E-74),  (B-75) , and  (B-76)  in  reference  3.  Thus 

/v 

Ya  and  by  noting  that  Q^(£»Y)  = NF  (y)  (s 


rt>  ± ( Y ) 


either  R or  Q we  obtain 

— n / m , a — n , m , o 


JF)  where  is 
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, , • 


J 


Mi(Y)  - (Ya) 


Pi(Y)  Polynomial  of  n degree 


We  will  now  assume  that  P^(Y)  has  simple  zeroes  (true  for  the 
sphere)  and  set 


N(i)  A . . 


i(y)  y 

k Y ^ 


From  equation  (93)  one  obtains 


b.  (yh  ^)  s • J . 

A.  . = -1  — 13  r "I:  e-^i(Yij) 


[dXi/dyl 
L J vzs'i 


and  equation  (89a)  can  be  rewritten  as 


«(i) 


J(r,Y) 


= V V e4»i(Y)-<fri(Yij)  

z-f  Z-f  e rdX./dYl 


i j = 1 


?i(^ij)  ^i^^ij) 

~7dyf  Y - Y±j 

jY=Yij 


For  the  sphere  <j>(Y)  = aY  and  consequently  equation  (95)  shows 
that  an  appropriate  rearrangement  of  the  eigenmode  solution  (89a) 
inexorably  leads  to  an  SEM  expansion  with  Class  1 coupling 
coefficients.  In  order  to  determine  the  form  of  the  entire  function 
(y)  for  a general  body  we  consider  the  inverse  Laplace  transform 
and  we  close  in  the  left  half-plane  for  times  t > t and  in  the 
right  half-plane  for  t < t where  tQ  is  the  instant  at  which  the 
incident  wavefront  first  hits  the  scattering  body.  All  the  poles 
are  located  in  the  left  half-plane  and  for  t < t we  should  obtain 
J(r,t)  = 0,  i.e.,  the  integral  on  the  large  semicircle,  Rey>o , should 
approach  zero  as  its  radius  approaches  infinity.  The  integrand 
of  interest  has  the  form 


,'nf* ' ..  ■•i.c-  > ■ v- 


AO 


( 

i 


1 

[ 


(96) 


(97) 

for  all  y along  the  semicircle.  Let  us  now  define  'F^(y)  by  the 
equation 

^(Y)  = (Y)  -yet  + yet 

e = t - t 

o 

For  t < t condition  (97)  can  be  rewritten 

Re  (y)  < | e | cRe  y . (98) 

Recalling  that  Re  y < | Y I we  can  cast  inequality  (98)  into  the 
form 


4^  (Y)  +Yct 


i (Y)  = 


Y - Y 


with  Re  Y^j  < 0.  Consequently  we  require  that 


Rej4K(Y)  + Yctj  < 0 


'J',  (Y) 

Re  cjTr < |e|  (99) 

with  e being  an  arbitrarily  small  finite  value.  Recalling  that 
4k(Y)  is  an  entire  function  one  would  be  tempted  to  expand  it 
in  a McLaurin  series  based  on  the  argument  that  equation  (99)  is 
only  required  to  be  true  on  a finite  semicircle  at  this  stage  of 
the  analysis  and  argue  that  equation  (99)  could  only  be  satisfied 
if  4>^(Y)  = constant  with  all  remaining  coefficients  in  the  series 
necessarily  equal  to  zero.  We  have  not  satisfied  ourselves  with 
this  argument;  however,  we  are  willing  to  conjecture  that  4>.(y) 

= constant  in  view  of  the  fact  that  any  polynomial  representation 
for  <|>^(y)  can  be  shown  to  reduce  to  a constant.  Making  this 
conjecture,  it  follows  that 
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(100) 


F=T 


i 
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( Y)  = -yet  + constant 

Such  a choice  of  course  allows  us  to  close  in  the  left-hand  plane 
for  t > t since  the  requirement 

Re  i^(Y)  < -ecRe  y = ec|Re  (Re  y < 0,  e > 0) 

is  then  satisfied  for  large  |y|  with  Re  y > 0. 

Assuming  that  4>^(y)  is  given  by  equation  (100)  then  <J>^(y)  - 
^i^Yij^  = ^Yij  ” Y ^ to  and  equation  (95)  is  an  SEM  solution  with  Class  1 
coupling  coefficients.  Rewriting  equation  (95)  after  making 
this  substitution  and  using  equation  (89),  one  obtains 


cii(^ 

J(r,Y)  = E(y)  £ 22 

i j-1  13 


(101a) 


where 


E (y)  = e~YCto 


(101b) 


cij(£) 


eYijcto  bi(Yij)  s ‘ di(^Yij) 


[dX./dyl 

L J Y=y 


1 D 


(101c) 


Combining  the  two  indices  (i  and  j)  into  the  index  "n,"  we  can 
write  equation  (101)  as 


V-'  c 

s • J(r,y)  = E (y)  ^ 


(102) 


n=l 


n 


and  this  form  is  convenient  for  a discussion  of  the  relationship 
between  meromorphic  functions  and  pole  series  expansions.  A 
consequence  of  the  Mittag-Lef fler  theorem  is  that  the  most  general 
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infinite  pole  series  representation  of  a meromorphic  function 
m(y)  is 


l(lf)  - E [p„(y^M  ' 5nH  + H<*> 


(103) 


where  and  Pn  are  finite  degree  polynomials  not  necessarily  of 

degree  n and  H(y)  is  an  entire  function.  There  exist  meromorphic 

functions  for  which  Pn  (l/ty  - Yn))  = Cn/ty  - Yn)  and  P (y)  and  H(y) 

are  identically  zero.  We  also  note  that  multiplying  the  meromorphic 

function  m(y)  by  an  entire  function  E(y),  would  still  result 

in  a meromorphic  function.  For  the  case  where  y and  C have  the 

n n 

appropriate  n dependence  to  eliminate  the  need  for  P^  and  H,  we 
would  obtain  a representation  for  the  meromorphic  function 
M(y)  = E(y)  m(y)  as  follows 


mm  = e(y) 


y - y. 


(104) 


which  is  consistent  with  both  equation  (102)  and  the  Mittag- 
Leffler  theorem.  This  is  also  consistent  with  our  assumed  initial 
representation  shown  in  equation  (89) . 
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SECTION  VIII 


NUMERICAL  DETERMINATION  OF  POLE  LOCATIONS  AND  EIGENFUNCTIONS 

FOR  THE  SPHERE 

In  this  section  we  describe  the  procedure  that  was  used  to 
approximately  find  the  pole  locations  for  the  sphere,  i.e.,  those 
values  of  y for  which  L(y)  - (1/2)1  has  a zero  eigenvalue.  To  do 
this  we  approximate  the  integral  operator  L(y)  -(1/2)1  by  a matrix, 
M(y),  according  to  the  algorithm  described  in  reference  7,  and 
try  to  find  those  values  of  y for  which  M(y)  has  a zero  eigenvalue. 
Since  the  matrix  is  of  finite  order  its  determinant,  D(y),  equals 
zero  if  and  only  if  at  least  one  of  the  eigenvalues  of  M(y)  is 
zero.  This  observation  permits  us  to  search  for  the  zeroes  of  D ( y ) 
rather  than  the  more  difficult  task  of  searching  for  the  zeroes 
of  the  eigenvalues. 

Our  method  of  calculation  of  the  elements  of  the  matrix  M(y) 
guarantees  that  D(y)  is  an  analytic  function  of  y.  Assuming  that 
the  errors  due  to  the  finite  word  size  of  the  computer  do  not 
interfere  with  certain  properties  of  analytic  functions  we  may 
apply  the  method  of  Singaraju,  Giri,  and  Baum  (ref.  8)  for 
determining  the  zeroes  of  an  analytic  function  to  find  those 
values  of  y for  which  D(y)  is  zero.  A brief  description  of  the 
method  follows. 

If  both  f(z)  and  g(z)  are  analytic  functions  of  z in  a 
simply  connected  open  domain,  ft,  then  Caughy's  theorem  states 
that  if  C is  the  border  of  a rectangular  region  contained  in  ft 


Sancer,  M.  I.,  S.  Siegel  and  A.  D.  Varvatsis,  Foundation 
of  the  Magnetic  Field  Integral  Equation  Code  for  the 
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Laboratory,  1976. 
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(105) 


' E 9<zi' 


where  the  are  the  N zeroes  of  f(z),  counted  according  to 
multiplicity , in  the  interior  of  C,  provided  that  f has  no  zeroes 
on  C.  In  particular  we  have  that 


N - 2ir/ ti 


(106) 


E <Vk  - 2ir/ 


k f • (z) 


k = 1,2, 


(107) 


and  from  these  equations  we  can  determine  the  z^. 

To  avoid  the  numerical  difficulties  associated  with  computing 
f'(z),  Singaraju  et  al.  (ref.  8 ) integrate  equation  (105)  by 
parts  to  obtain 


E gui' = 


N g(zinit)  + g'(z)  adj  log  <f (z) ) dz  (108) 

c 


where  adj  log  (f(z))  is  that  branch  of  the  locally  defined  log 
function  whose  only  discontinuity  along  the  contour  occurs  at 
zinit  where  undergoes  a jump  of  2ffi  N.  The  problem  of 


Singaraju,  B.  K.,  D.  V.  Giri  and  C.  E.  Baum,  Further 
Developments  in  the  Application  of  Contour  Integration  to 
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Relevant  Computer  Programs,  Mathematics  Note  42,  Air  Force 
Weapons  Laboratory,  1976. 
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calculating  f ' (z)  has  been  replaced  by  the  simpler  problem  of 
finding  the  phase  of  f(z)  and  guaranteeing  that  the  phase  change 
from  point  to  point  is  less  than  tt  radians.  Using  equation  (108), 
equations  (107)  become 


N 


5 E 

i=l 


z . = 


k 

!init 


>TT1  J 


Zk  ^ adj  log(f(z))  dz  k=l,2. 


. . ,N 
(109 


For  further  details  we  refer  the  reader  to  reference  7. 

In  their  paper,  Singaraju  et  al.  provide  a computer  program 
to  determine  the  z^  if  a contour  contains  up  to  three  zeroes, 
however,  the  sphere  problem  has  a high  degeneracy  and  therefore 
the  numerical  approximation  has  the  zeroes  of  D(y)  clustering 
about  their  true  locations.  We  therefore  found  it  necessary  to 
extend  their  procedure  to  arbitrary  N.  To  accomplish  this  we 

i.  u. 

transformed  the  system  of  equations  (109)  to  the  equivalent  N 
degree  polynomial  equation 


N 


k=0 


(110- 


where  are  defined  by  the  recursive  formula 


Bo  " 1 


(110) 


kBk  + 


T:  Aj  Bk_j  = 0 k = 1,  ...,  N 

j = l 


Equation  (110b)  was  found  in  reference  9 as  suggested  by  Dr.  Giri 


i 


i 

* 

I 

'* 
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9.  Krylov,  V.  I.,  Approximate  Calculation  of  Integrals, 


(Priblizhennoe  Vychislenie  Inteqralov) , Moskva,  Gos.  Izd-Vo 


We  have  chosen  to  use  the  subroutine  package  ZRPCC  of  the 
AFWL  Scientific  Program  Library  (ASPLIB)  to  find  the  zeroes  of  this 
polynomial.  This  routine  forms  the  companion  matrix  whose  charac- 
teristic equation  is  identical  to  the  required  polynomial  and  then 
returns  the  eigenvalues  of  that  matrix. 

For  each  y that  estimates  a zero  of  D(y)  and  therefore  a pole 
location  of  the  sphere's  response  function  we  compute  a new  M(y). 

We  "feed"  this  matrix  to  another  ASPLIB  subroutine  package,  CIVAA, 
which  is  documented  as  CG  in  the  Eispack  Guide  (ref.  10).  This 
package  computes  the  eigenvalues  and  eigenvectors  of  M(y).  We 
expect  to  find,  and  in  practice  have  found  that  all  of  the  eigen- 
values occur  in  pairs  which  sum  to  one  (see  equation  (24  )). 

We  search  the  list  of  eigenvalues  for  those  pairs  near  (0,1) 
and  identify  the  corresponding  eigenvectors  as  the  natural  mode 
vector  J and  its  dual  vector  J.  In  theory,  because  we  are  solving 
the  problem  for  a sphere,  J and  J are  related  by  the  formula 

a /-w»  ^0  i 0 * 

J(r)  = n(r)  x J(r)  e where  the  e is  introduced  to  account 
for  the  arbitrariness  in  the  phase  of  an  eigenvalue;  our  numerical 
experiments  showed  that  this  relationship  held  to  within  the  errors 
of  our  approximation. 

Since  the  sphere  problem  has  a high  degeneracy  we  cannot  expect 
our  eigenfunctions  to  agree  with  those  published  in  other  sources 
but  we  can  expect  ours  to  be  linear  combinations  of  theirs.  As 
an  example,  for  the  n = 1,  q = 1 mode  we  found  three  eigenmodes 
J^,  J_2>  £3  which  can  be  related  to  those  found  in  reference  2. 
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APPENDIX  A 


PROOF  OF  EQUATION  (14) 


•j*  *j* 

In  this  appendix  we  show  that  Q = Q*  and  M 


= -M , i . e . , 


equation  ( 14  ) is  true.  Recall  that 


QJ  (r)  = n (r)  x (n  ( 


r)  x j f (R) (r  - r')  * J(r 


')  dS') 


/ 


= -P (r)  • / f(R)(r  - r’ ) * J(r')  dS’ 

'S 


where 


P (r)  = I - n(r)  n(r)  = -n(r)  x ( (n  (r)  x i] 

f (R)  = -<1  + YR)(^),  R = |r  - r'| 

Notice  that  for  a tangential  vector  a,  P • a = a.  Consider  now 
the  inner  product  involving  the  tangential  vector  functions  Mr) 
and  41  (r)  : 


(*,QY)  s f **(r' 


) • (QH'(r))  dS 


= - j" dS | <5 * ( r ) • P (r)  • J jf (R) (r  - r ’ ) x Y(r ’ )j  dS'j 


-// 


dS  dS’  f (R) 


j{^*(r)  x (r  - r')J  • Mr ' )j 


If  we  interchange  r and  r'  the  "double"  integral  does  not  change 
value  and 

*9 


(* 


*,qi)  = -JJ  dS'  dS  f(R)j^£*(r')  x (r1  - r)J  • ¥(r)  > 
S'S  ^ * 


= - j J dS  ' dS  [f  (R)  (r  - r ' ) X * * (r  ' )J  • ¥ (r) 
S:S 


= - / dS  ¥(r)  • P (r)  • i ["f*(R)(r  - r’)  x £(r')] 

Q *c  » 


= JdS  j-P(r)  • J f*  (R)  (r  - r ' ) x£(r')j 


dS' 


/< 


= (Q*  <P  (r)  ) • Y(r)  dS 


(Q*  *,v) 


Recalling  the  defining  relationship  for  we  conclude  that 
Q+  = Q*. 

Next  we  show  that  = -M. 


(£,MY)  5 j 


dS  <I>*  • (n 


x v • - J 


dS  (n  x $*)  . v 


■/ 


dS  (-n  x $)  • Y = (-M£,£) 


and  this 


proves  that  M 
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APPENDIX  B 


PROOF  THAT  Nj_  dX±/dy  = Jd • |n  x jJ 


In  this  appendix  we  show  the  validity  of  the  formulas  employed 
in  section  V for  the  coupling  coefficients  by  demonstrating  that 


'i  ^T=/ds[ii  • (”  * If  ii)] 


(B-l ) 


where  we  have  ignored  degeneracy  since  the  proof  is  identical  for 
the  degenerate  case.  We  start  with  equation 


- h h 


and  differentiate  both  sides  w.r.t.  y 


r\Cf  • 

j.  + 3?-—  = j • + X . . 

9y  —l  9y  dy  — l l 9y 


We  operate  on  both  sides  by  M and  consider  the  inner  product  of 
the  resulting  equation  with  FV  defined  in  equation  (9). 


(li-  « |f  -i)  * (Fi-  = 1FT  + xi(Fi'  M Tr) 


(B-2) 


By  definition  (FV,  M.f/ 3j^/9y)  = ((M'/O^F.,  9J^/9y)  and  from  equation 

(18)  = Ft.  We  can  now  rewrite  equation  (B-2)  as 
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